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Transport of molecules across membrane channels is investigated theoretically us- 
ing exactly solvable discrete stochastic site-binding models. It is shown that the in- 
teraction potential between molecules and the channel has a strong effect on translo- 
cation dynamics. The presence of attractive binding sites in the pore accelerates the 
particle current for small concentrations outside of the membrane, while for large 
concentrations, surprisingly, repulsive binding sites produce the most optimal trans- 
port. In addition, asymmetry of the interaction potential also strongly influences 
the channel transport. The mechanism of these phenomena are discussed using the 
details of particle dynamics at the binding sites. 

Membrane proteins support and regulate fluxes of ions and molecules that are critical 
for cell functioning 1]. Molecular transport across membrane pores is characterized by high 
efficiency, selectivity and robustness in the response to fluctuations in the cellular environ- 
ment, however precise mechanisms of these complex processes are still not well understood 
Membrane channel proteins, that utilize the energy of adenosine triphosphate 
(ATP) hydrolysis to move ions and small molecules against external free energy gradients, 
are known as active transporters. There are experimental and theoretical arguments that 
suggest that the high selectivity in these proteins is reached via specific interactions at the 
narrowest part of the pore although some recent experiments indicate that non-specific 
interactions with the whole channel are also important Q|. It was assumed earlier that the 
membrane proteins with large water-filled pores move molecular solutes in a passive trans- 
port mode by utilizing a simple diffusion, and that these large channels have relatively low 
efficiency and selectivity. However, recent experiments suggest that permeating molecules 
interact strongly with large membrane pores leading to a very efficient and highly-selective 
transport BflflflQ. 

In order to understand the facilitated transport phenomena in large membrane pores sev- 



eral theoretical approaches have been presented Q, [ill E3, Q Q 

A continuous model that 

describes the motion of a single solute molecule in the channel as one-dimensional diffusion 
a,on g the potential of ,„ea n forces with position-dependent diffusion constant QQQ 
has investigated the most efficient permeation dynamics, and it was shown that there is an 
optimum attraction between the channel and the translocating molecule that creates the 
maximal flux across the membrane. However, a uniform potential along the entire channel 
with the attraction magnitude of 6-8 /c^T has been assumed in the calculations that repro- 
duce the experimentally measured currents. It is known that the structure of membrane 
channel proteins is very complex Q], and it is reasonable to suggest that the corresponding 
free-energy potential of interaction has a very rough landscape. Recent molecular dynamics 
simulations of glycerol translocation through aquaglyceroporin GlpF [l^ calculate a po- 
tential of mean forces that shows several relatively weak (2-6 ksT) but strongly localized 
narrow wells. In a different approach, a macroscopic version of Fick's law has been used to 
analyze the molecular transport through the membrane channels and it was concluded 
that in the idealized model any interaction would lead to the amplification of molecular flow 
through the channel. However, this conclusion is rather unphysical since the binding site 
with infinite attraction would block the molecular traffic through the channel, contradicting 
to this theoretical prediction. 

Potentials of interactions between the channels and permeating molecules, determined 
in experts and from confer sanations Q |J are g ene rally as y „rmetrrc with 
multiple weakly attractive and repulsive binding sites. Although the coupling between 
non-equilibrium fluctuations and asymmetric potential of interactions in channels has been 
discussed recently \^\ . a full theoretical description of the effect of attractions and repul- 
sions, and the asymmetry of the interaction potential on the molecular currents through 
the biological pores is not available. In this paper I present a theoretical analysis of the 
channel-facilitated transport using simple discrete stochastic models with multiple binding 
sites that allow to calculate exactly the dynamic properties of the system. It is found that 
depending on the concentrations outside of the membrane attractive or repulsive binding 
sites might increase molecular fluxes, and the translocation dynamics is strongly influenced 
by the asymmetry of the potential. 

I consider the transport of particles across a cylindrical membrane channel as shown in 
Fig. 1. The molecule enters the pore from the left (right) chamber that has a concentration 



3 



ci (c 2 ) with the rate u = k on ci (w = fc on c 2 ); and it leaves the channel to the left (right) with 
the rate W\ = k ff (u^ = k fj). It is assumed that particles are not interacting with each 
other, and there are N generally non-equivalent binding sites in the channel. The particle 
at the site j (j = 1, 2, • • • , N) jumps to the right (left) with the rate Uj (wj). Defining Pj(t) 
as a probability of finding the molecule at the binding site j at time t, the translocation 
dynamics can be described by a set of master equations, 

dPj(t) 



dt 



Uj^Pj^t) + W j+1 P j+1 (t) - (Uj + Wj)Pj(t), 



(1) 



where j = 1, 2, • ■ • , N an d Pn (t) = P/v + i(t) = 1 — ^ x Pj(t) is the probability of finding the 
channel empty at time t 13|. The discrete stochastic model with N binding sites can be 
solved exactly by mapping it into a single-particle hopping model along the periodic infinite 



one-dimensional lattice with N + 1 states per period 1£, [3]. This model has been used 
successfully to describe the dynamics of motor proteins Ji| . The mapping can be easily seen 
by considering multiple identical channels from Fig. 1 arranged sequentially and keeping 
the concentration gradient across each period as Ac = C\ — c 2 . Then dynamic properties of 
the channel-facilitated transport model can be calculated exactly. Specifically, the particle 
current J for the system shown in Fig. 1 with N identical sites without interactions (uj = 
Wj + i = a for j = 1, 2, • ■ ■ , N — 1) is given by 13 
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(2) 



To investigate the effect of interactions in the membrane transport the simplest model 
with N = 1 binding site is analyzed. Assume that energy of the binding site is equal to — e, 
i.e., e > corresponds to the attractive site, while negative e describes the repulsive site. 
The transition rates are related via detailed balance conditions, 



uq(e) u (e = 0) Ui(e) U\(e = 0) 



w±(e) Wi(e = 0) w (e) w (e = 0) 
and they can be written in the following form 



1/x), with x = exp (e/fc B T), 



(3) 



uo(e) = uqx 01 , wi(e) = w\x dl 1 , U\(e) = uqx 02 1 , Wq(s) = wqx 9 ' 2 , 



(4) 



where interaction-distribution coefficients d\ (with < 9i < 1 for % = 1, 2) describe how the 
potential changes the transitions rates. For simplicity, it is assumed that 9\ = 9 2 = 9. Then 
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the particle current across the membrane channel is equal to 

{UqUi - W Wi)x e k n(Cl-C 2 )x 6 



J 



(u + W )x + (Ui + W\) 2 -f fc °"( c i+ c 2) 



(5) 



X 



Define Jo as the molecular flux in the system without interactions (e = 0), then the dimen 
sionless ratio of currents is given by 

J _ [kon(ci + c 2 ) + 2k off ] x 6 



(6) 



J 2k off + k on {c x + c 2 )x 

The molecular flux across the membrane depends strongly on the interaction strength at 
the binding site, as shown in Fig. 2 for different concentrations. For strong attractions 
and repulsions the current decrease, while for intermediate interactions the molecular flow 
increases significantly, reaching the maximum value at e*, 

9 2k off 



e* = k B T\n 



(7) 



[1 - 9) k on (d + c 2 

At this interaction the molecular transport across the membrane is the most optimal with 



the most efficient relative current 



T ) =(1-*) 
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off 



1 - 9) k on (d + c 2 ) 



(8) 



It can be shown that the largest increase in the molecular flux can be achieved when 9 — > 1, 
producing 



l+ ; ; • (9) 

kon{Cl + C 2 )_ 

The effect of current increase can be estimated explicitly for membrane channel proteins by 
utilizing the transition rates k ff ~ 500 s _1 and k on ~ 15 yuM _1 s _1 from the experiments on 
maltodextrin translocation through maltoporin channels [9|. For concentrations of few /iM 
Eq. © predicts r^IOO times increase in the molecular fluxes. 

The most optimal interaction e* depends on the molecular concentrations outside of 
the membrane pore as illustrated in Fig. 3. Our analysis shows that the presence of the 
attractive site leads to molecular flux increase for small concentrations C\, while for large 
concentrations the presence of the repulsive site increases the particle current, and there 
is a critical concentration c* (for every fixed value of the concentration c 2 ) that separates 
two regimes. The fact that the repulsive binding site leads to the molecular flux increase 
seems, at first, surprising, and it was not predicted by previous theoretical approaches 



jlQ.lllUl2lll3l.ll4l]. However these observations can be explained using following arguments. 
Assume that the concentration to the right of the membrane is zero (see Fig. 1), i.e., = 0, 
although our arguments can be easily generalized. It can be shown that the current across 
the membrane can be viewed as a ratio between the effective probability to translocate the 
pore over the mean residence time that particle spends in the channel |20J], J = U/t. In this 
case II = 1, and the current is inversely proportional to the translocation time r [soj ] . 

1 1 wi 2 x 1 - 8 1/2 x \ 

T = 0+ + 6=1 6+1 = -fl 1 + 1 ' 10 

U X y UiX"- 1 UqUiX K nC\X ti k of f X° \k on Ci k of f J 

This expression can be understood as a sum of two contributions, namely, the effective 
times to enter the binding site and to leave it. The conditions of the most optimal transport 
correspond to the situation when these two terms are approximately equal. Then the increase 
in c\ lowers the value of the most optimal interaction, and for large concentrations the 
repulsive binding site provides the most efficient translocation. It is reasonable to suggest 
that in the general case of interaction potential with N non-equivalent binding sites the most 
optimal transport is achieved when the effective times to enter and to leave the strongest 
attractive or repulsive sites balance each other. 

The potential of interaction between the solute and the membrane channel is asymmetric 
jlfl llfi| . and to study this property of the molecular transport a discrete stochastic model 
with N = 2 binding sites is used. The asymmetry is introduced by assuming that energies 
of two consecutive binding sites are — e and or and — e, respectively. If the interaction is 
placed on the first binding site, the particle current across the membrane is given by 

k Jci - c 2 ) 



I I fcongl I x -6 M _|_ k °ff _|_ fcon c 2 I k on C2 \ I rjA-9 I k on C2 _i_ kon£l \ _|_ rv. k on C\ 
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while for the case when the interaction is on the second binding site one can obtain 



-y _i_ k on C2 _i_ M _)_ k °ff _)_ k „ci , k on ci \ , x l-8 ( k on c\ , k on C2 \ _i_ ^koncg 

k ff ^ { ^ a ^ koff a ) ^ \ k oft ^ Ot J k ff 



;i2) 



where the same interaction-distribution factors 9 are assumed in both cases. The ratio of 
two currents, plotted in Fig. 4, deviates from unity for all interactions except e = 0. For 
very strong attractions (e — > 00, x ^> 1) one can show that 

~ — <1, for c 2 < ci, (13) 
J2 Ci 
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while for strong repulsions (e — > — oo, x — > 0) 



, ^ _|_ fcongl 

fco// " >1, forc 2 < Cl . (14) 
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Generally, for e < we have J x / J 2 > 1, while for attractive interactions the ratio of currents 
is always smaller than one. Thus putting the binding site at different positions along the 
channel changes the molecular flux across the membrane. This surprising observation can 
be explained by looking at dynamics of the particle translocation across the channel. First, 
consider the repulsive interaction at the first binding site. After the particle passes the 
binding site it has a low probability to come back from the second site since the barrier 
is high. As a result, the overall translocation time is low and the translocation current is 
large. However, when the repulsive interaction is at the second binding site, the situation 
is different. After reaching the binding site the particle has a high probability to return to 
the first site, and many attempts to cross the second binding site will be made before the 
successful passing the channel. As a result, the overall translocation time is high, leading to 
small molecular fluxes across the pore. 

Our theoretical approach is based on the discrete-state stochastic description of molecular 
fluxes through membrane channels, although the majority of other theoretical approaches 
have utilized continuum models [ill 13, (liil Q|- There are several advantages of using 
discrete models of membrane transport, namely, the ability to describe better the nature 
of interactions between the membrane proteins and translocating molecules, biochemical 
complexity of molecular transport and existence of exact and explicit results for dynamic 
properties Q]. In addition, it was shown recently that continuum and discrete models of 
membrane transport produce equivalent results 13j | . 

In conclusion, discrete-state stochastic models of translocation of molecules across mem- 
brane pores, that allow to calculate explicitly the dynamic properties of the system, are 
presented. The conditions for the most optimal channel transport are discussed. It is shown 
that the strength of interactions at the binding sites strongly influences the translocation 
dynamics. For small concentrations outside of the membrane the attractive sites yield the 
largest particle current, and the repulsive binding sites produce the most efficient molecular 
transport for large external concentrations. It is argued that the most optimal conditions for 
the membrane transport are achieved when the times to enter attractive or repulsive sites 
are balanced by the corresponding times to leave these positions. To best of our knowledge, 
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this work is the first that suggests that repulsive interactions might be favorable for translo- 
cations across the channels. Our theoretical method is also used to investigate the effect of 
asymmetry on the membrane transport by putting a strong interaction site at different posi- 
tions in the channel. The asymmetry gives different free-energy landscapes, thus producing 
different molecular currents. This suggests that the asymmetry in the interaction potential, 
strongly effects the overall membrane transport, even without coupling to non-equilibrium 
fluctuations [l7J. The presented theoretical analysis supports the idea that interaction be- 
tween the molecules and the entire channel is important for the selectivity and efficiency of 
the membrane transport J|. 

Our calculations indicate that the mechanism of high selectivity and efficiency of mem- 
brane channels is due to the interaction potential between molecules and the channel, and it 
is reasonable to suggest that the evolution tuned this potential to create the most optimal 
molecular transport It is important to note, however, that our model is rather oversim- 
plified and many important factors of translocations dynamics, such as interactions between 
molecules and three-dimensional nature of the channels and corresponding interaction po- 
tentials, are neglected. At the same time, it is expected that the main physical principles of 
the translocation across membrane pores presented in this work are still generally valid. It 
will be important to investigate the validity of these predictions by analyzing more realistic 
models of membrane transport. 
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Figure Captions: 

Fig. 1. A general kinetic scheme for a stochastic model of membrane transport with TV 
binding sites. A membrane separates two chambers with concentrations c\ and ci- A 
particle can enter the channel from the left with the rate Uq = k on ci or from the right with 
the rate w = k on C2, and it leaves the pore with rates Un — W\ — fc //- At the site j the 
particle jump forward and backward with rates Uj and Wj, respectively. 

Fig. 2. Relative currents as a function of the interaction strength for the model with N = 1 
binding site for different concentrations. The transitions rates, k on = 15 /iM^'s -1 and 
k ff = 500 s" 1 are taken from Ref. For all calculations C2 = is assumed. Different 
curves correspond to calculations using Eq. © with a) C\ = 10 /iM and 9 = 0.5; b) c\ = 10 
/M and 6 = 0.9; c) d = 500 /iM and 6 = 0.5; and d) c x = 500 /iM and 9 = 0.9. 

Fig. 3. The most optimal interaction as a function of the external molecular concentration 
C\. For calculations C2 = and 9 = 0.5 are assumed. 

Fig. 4. The ratio of two currents (see text) as a function of the interaction strength for the 
model with N = 2 binding sites. The transitions rates, k on = 15 /iM _1 s _1 and k ff = 500 
s _1 are taken from Ref. [9j. For all calculations C2 = 0, 9 = 0.5 and a = k ff are assumed. 
Different curves correspond to different concentrations: a) c\ = 10 /iM; b) c\ = 100 /iM; 
and c) a = 300 /iM. 
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Figure 1: A. Kolomeisky, Physical Review Letters. 



11 




zlk B T 



Figure 2: A.B. Kolomeisky, Physical Review Letters. 
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Figure 3: A.B. Kolomeisky, Physical Review Letters. 
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Figure 4: A.B. Kolomeisky, Physical Review Letters. 



